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ABSTRACT 

We report the detection of radio QPOs with structure changes using the Very Long 
Basehne Array (VLB A) at 43 GHz. We found conspicuous patterned changes of the 
structure with P = 16.8 ± 1.4, 22.2 ± 1.4, 31.2 ± 1.5, 56.4 ± 6 min very roughly in a 
3 : 4 : 6 : 10 ratio. The first two periods show a rotating one-arm structure, while the 
P = 31.4 min shows a rotating 3-arm structure, as if viewed edge-on. At the central 
50//as the P — 56.4 min period shows a double amplitude variation of those in its 
surroundings. Spatial distributions of the oscillation periods suggest that the disk of 
SgrA* is roughly edge-on, rotating around an axis with PA = —10°. Presumably, the 
observed VLBI images of SgrA* at 43 GHz retain several features of the black hole 
accretion disk of SgrA* in spite of being obscured and broadened by scattering of 
surrounding plasma. 
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1 INTRODUCTION 

The existence of black holes has been definitely established 
(Miyoshi et al. 1995; and Harms et al. 1994) while zooming- 
in the relativistic region is still in difficulty (Tanaka et al. 
1995) though promising in near future (Falcke et al. 2000, 
Miyoshi et al. 2004, 2007, Doeleman et al. 2008, Fabian 
2009). Sagittarius A* (SgrA*), the most convincing massive 
black hole (Ghez et al. 2001; Schodel et al. 2002, and Shen et 
al. 2005) at the Galactic center, shows short time flares with 
quasi-periodic oscillations (QPO) with P = 17, 22, & 33 min 
in near-infrared and X-ray regions originated from near the 
central black hole (Baganoff et al. 2000; Genzel et al. 2003; 
Aschenbach et al. 04; Eckart et al. 2006; Belanger et al. 2006; 
Hamaus et al. 2009) 

However the infrared and X-ray observations of SgrA* are 
limited to one or two hour's durations of the flare events. 
Due to the short time durations, the reliability of the QPO 
periodicity is still open to argument. Also the dependence 
of the detection upon observing instruments is another ar- 
gument. (Meyer et al. 2008) 

In contrast, because the SgrA*is bright all time at ra- 
dio wavelength, the continuous observations about 7 hours 
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can be performed even from the northern hemisphere. The 
longer durations of radio observations of SgrA* give us a 
chance to enhance the reliability of the QPO periodicity if 
the QPOs occur at radio wavelength. Further, if we can ob- 
serve the QPO periodicity of SgTA*not only in intensity but 
also in structure change, the reliability of QPO detection 
will be conclusively enhanced. In practice, however it has 
been recognized that the calibration of the VLBI data on 
SgrA* is quite difficult because of the atmospheric variations 
at the low elevation observations (Bower et al. 1999). Then 
the VLBI synthesis imaging of SgrA*is recognized as quite 
difficult, therefore the estimations of the size and shape of 
SgrA*were performed through the closure analysis that are 
free from instrumental and atmospheric errors (Doeleman 
et al. 2001, Bowers et al. 2004, Shen et al. 2005) As the 
closure quantities are not so robust against thermal noise, 
the use is limited to high SNR visibility data from shorter 
baselines (less than 2000 km). Though the corresponding 
spatial resolutions or synthesized beam sizes are somewhat 
larger than the size of SgrA* at 43GHz, the closure analysis 
can be performed successfully at the observing frequency. 
Recently Miyoshi et al. (submitted to PAS J), by devising 
a simple method at calibration, effectively used the data 
from longer baselines, and achieved reliable synthesis imag- 
ings of SgrA* at 43GHz, which are at least consistent to the 
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previous closure estimations. Using the VLBA, we inves- 
tigated the short time change with another new method, 
slit-modulation imaging method (Miyoshi 2008), and found 
four of periods with conspicuous patterned changes of the 
structure. 



2 OBSERVATION AND DATA CALIBRATION 

Our VLBA observations at 43 GHz were performed with 512 
Mbit-per-second recording rate (16 MHz x 8 IF channels, 2- 
bit sampling) from 9 : 30 to 16 : 30(UT) on 2004 March 8"*, 
37 ~ 44 hours just after a millimeter wave flare of SgrA* de- 
tected with the Nobeyama Millimeter Array (Miyazaki et al. 
2006). 

In order to make efficient use of as much data as pos- 
sible, we adopted a reduction strategy of keeping all data 
points till the final calibration stage where the CALIB (AIPS 
task of self-calibration) is performed. For data points where 
no solution was obtained at FRING in AIPS, wc assigned 
a value interpolated from adjacent good solutions to the 
data points. Such keeping strategy in reduction is often per- 
formed for weaker maser data. Wc did the standard am- 
plitude calibrations using antenna gains in GC table and 
system temperatures in TY table with opacity corrections 
in AIPS (NRAO). Further at the last stage of the data cal- 
ibrations we applied solutions of amplitude and phase from 
the task CALIB with an image model of the average size and 
shape of SgrA* at 43 GHz (single elliptical Gaussian with 
full width at half maximum of 712/ias x 407^ias, PA = 79.8°, 
Bower et al. 2004) with 1 Jy in flux density. The application 
of the CALIB solutions to the data substantially reduced 
the effect of atmospheric and instrumental variations on the 
observed visibilities (Miyoshi et al. submitted to PAS J), and 
effectively made visibilities of longer baselines available for 
imaging synthesis. Figure 1 shows that the u-v coverage of 
the solutions extends up to 5800 km (8.3x 10* A) in projected 
baseline length. This will give us a spatial resolution of three 
times better than those of the previous studies (Bower et al. 
1998, 2004). We show the amplitude and phase variations of 
the calibrated visibilities of respective baselines and closure 
phases for all triangles in Figures 2, 3, & 4 and Figures 5, 
6, & 7 respectively. 

2.1 The u-v coverage of the calibrated data 

We got calibration solutions from 43 baselines out of all the 
45 baselines, though not during the whole observing time. 
While we found no solution from SC-MK and HN-MK base- 
lines. As for the baselines connected to MK, we found effec- 
tive visibilities when the projected baseline length becomes 
shorter than 5800 km during the last one-half hours of the 
observing session. Except SC-MK and HN-MK, we found so- 
lutions from the 7 baschnes (BR-MK, FD-MK, KP-MK, LA- 
MK, NL-MK, OV-MK, and PT-MK). As for the baselines 
connected to SC, we found solutions from the 8 baselines 
(BR-SC, FD-SC, HN-SC, KP-SC, LA-SC, NL-SC, OV-SC, 
and PT-SC). We found none of solution from the SC-MK 
baseline. As for the baseline connected to HN, except the 
baseline between MK we found solutions from the 8 base- 
lines (BR-HN, FD-HN, HN-KP, HN-LA, HN-NL, HN-OV, 
HN-PT, and HN-SC). As for the baseUnes connected to BR, 



we found solutions from all the baselines. Though the BR 
station locates at high latitude, we often find good visibili- 
ties of SgrA* because the projected baseline length becomes 
shorter for the SgrA*observations. The all of resultant u- 
V coverage distribute up to 8.3 x 10*A (5800 km) in radius 
(Figure 1). 

2.2 Visibility variation in each baseline 

Figures 2, 3, & 4 show the time variations of the amplitude 
and phase of each baseline (each point is 5 min integration) . 
The lines in every plot show the time variations of the am- 
plitude and phase calculated from the obtained image from 
the whole time integration (Figure 8 a). From the compar- 
ison between the real visibilities and calculated ones from 
the image, we can classify them into 6 categories. 

(i) The phase is constantly zero and the amplitude has 
no variation with time. Namely the source SgrA*is observed 
as an unresolved point source for the corresponding fringe 
spacing of the baselines: 2 baselines, KP-PT and LA-PT 
show this feature. The maximum projected baseline lengths 
are about 400 km (KP-PT) and 210 km (LA-PT). The cor- 
responding fringe spacings are 3.5 mas and 6.7 mas respec- 
tively. 

(ii) The calculated phases and amplitudes from the ob- 
tained image, and the real ones coincide each other well. 
The phases are constantly zero, while the amplitudes vary 
with projected baseline length. The observed source SgrA* is 
partially resolved by the spatial resolutions of the baselines. 
In other words, SgrA*was observed to be not infinitely small 
but have a size: 19 baselines, namely BR-FD, BR-KP, BR- 
LA, BR-OV, BR-PT, FD-KP, FD-LA, FD-NL, FD-OV, FD- 
PT, HN-NL, HN-SC, KP-LA, KP-NL, KP-OV, LA-NL, LA- 
OV, NL-PT, & OV-PT, belong to this category. The aver- 
age of these maximum projected baseline lengths is about 
1400 km (1 mas in fringe spacing: from here in bracket after 
projected baseline length, we note the corresponding fringe 
spacing), the shortest ones are about 560 km (2.5 mas) (FD- 
LA, FD-PT) and the longest one is about 2650 km (0.53 
mas) (HN-SC). 

(iii) The calculated phases from the obtained image, and 
the real ones coincide each other or show similarity. While 
the amplitudes from both do not show similarity. 3 baselines, 
BR-SC, BR-MK, PT-MK, belong to this category. The max- 
imum projected baseline lengths are 5740 km (0.24 mas), 
3850 km (0.36 mas), and 4760 km (0.29 mas) respectively. 

(iv) The calculated phases from the obtained image, and 
the real ones do not coincide each other. While the am- 
plitudes from both show similarity and match in some de- 
gree: 12 baselines, BR-HN, BR-NL, FD-HN, HN-KP, HN- 
LA, HN-OV, HN-PT, KP-SC, NL-OV, HN-SC, OV-SC, & 
PT-SC, belong to this category. The average of these max- 
imum projected baseline lengths is about 3710 km (0.38 
mas), the shortest one is about 2380 km (0.59 mas) (NL-OV) 
and the longest one is about 5530 km (0.25 mas) (OV-SC). 

(v) The calculated visibility and the real ones do not show 
matching both in phase and in amplitude: 7 baselines, FD- 
SC, FD-MK, KP-MK, LA-SC, LA-MK, NL-MK, & OV-MK, 
belong to this category. In general there are at least three 
possibilities to explain the no-matching. One is due to a low 
signal to noise ratio including the case of no signal. Another 
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Figure 1. The u-v coverage of all the visibility data that have calibration solutions by our reduction method. 



is due to the complex structure of the real source, the ob- 
tained rough image could not reproduce the real visibiUty 
variations. The other is the mixture of the two cases above. 
The average of these maximum projected basehne lengths is 
about 4650 km (0.30 mas), the shortest ones arc about 3920 
km (0.36 mas) (NL-OV) and the longest one is about 5670 
km (0.25 mas) (OV-SC). 

(vi) No solution is found in SC-MK and HN-MK base- 
lines, which are longer than 5800 km in projected baseline 
length. 

Assuming that matching between the observed and the 
calculated visibilities mean sufficient calibrations of the two 
stations of the baseline, we can guess whether the calibra- 
tions both of the stations are good or bad. (Here we think 



that mis-matching does not necessarily mean bad calibra- 
tions. If the observed visibility is from time variable source 
in structure, any single image cannot satisfy the whole time 
visibility variation. However, if the structure variation is not 
so large one, such a small variation is negligible and not de- 
tectable for shorter baselines.) By eye, we found that the 
following baselines seem in the matching: BR-FD, BR-KP, 
BR-LA, BR-OV, BR-PT, FD-KP, FD-LA, FD-NL, FD-OV, 
FD-PT, HN-NL, HN-SC, KP-NL, KP-OV, KP-PT, LA-NL, 
LA-OV, LA-PT, NL-PT, and OV-PT. If the assumption is 
correct, except the MK station, other 9 stations seem to be 
well-calibrated anyway. 
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Figure 2. The time variations of visibility data in respective baselines (1/3). 
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Figure 3. The time variations of visibility data in respective baselines (2/3). 
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Figure 4. The time variations of visibility data in respective baselines (3/3). 
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2.3 Closure phases 

Figures 5, 6, & 7 show the time variations of the closure 
phases of all triangles composed from three stations. Every 
point is from 15 min integration. We also calculated the 
closure phase variations from the images of the whole time 
integration as shown in Figure 8 (a). Closure phases show 
several features as below. 

(i) First feature is seen in the small triangles composed 
from inner stations of the VLBA, namely FD, LA, FT, KP, 
& OV. The projected extensions of the triangles range from 
700 km (2 mas) to 1600 km (0.88 mas), and the average 
is 1200 km (1.2 mas). Those 10 triangles show nearly zero 
constant closure phases, a few degree variations at most. 
This indicates that SgrA*is observed as a structure of point 
symmetric with the corresponding fringe spacings (0.88 — 
2mas) . 

(ii) Second feature is seen in the 10 triangles composed 
from BR and other two stations from the inner 5 stations, 
FD, LA, FT, KP, & OV. The projected extensions of the 
triangles range from 1500 km (0.93 mas) to 2200 km (0.64 
mas), and the average is 1900 km (0.74 mas). The closure 
phases show zero on average, but fluctuating ±40° at peak. 

(iii) Third feature is seen in the 10 triangles composed 
from NL and other two stations from the inner 5 stations, 
FD, LA, FT, KP, & OV. The projected extensions of the 
triangles range from 1800 km (0.78 mas) to 2800 km (0.5 
mas), and the average is 2400 km (0.58 mas). The closure 
phases show zero on average, but fluctuating by ±60° on 
average, ±100° at maximum case. 

(iv) In most of large triangles which have projected exten- 
sions up to 5800 km (0.24 mas), the closure phases distribute 
±180°. However, some of triangles including MK station 
show nearly constant and/or very small variations in closure 
phase though the duration is about 40 minutes. Two triangle 
of FD-NL-MK and L A-PT-MK show closure phase variation 
within 15°, triangle of KP-PT-MK shows that within 10°, 
and triangle of KP-LA-MK shows nearly constant closure 
phase within 8° variation. The projected extensions of these 
three cases range from 4800 km (0.29 mas) to 5200 km (0.27 
mas), 4900 km (0.29 mas) respectively. 

(v) In most of the large triangles, the real closure phase 
and the calculated phase from the obtained image do not 
match. While one of the large triangles, HN-NL-OV shows 
matching between them. The projected extension of the tri- 
angles is about 4200 km (0.33mas). 

(vi) The last feature is seen in the five of triangles BR- 
FD-KP, BR-FD-LA, BR-FD-NL, BR-FD-OV, BR-FD-PT, 
whose projected extensions range from 2100 (0.67 mas) km 
to 2500 km (0.56 mas), and the average is 2300 km (0.64 
mas). They show very similar closure phase variations each 
other. 

2.4 The calibrations of MK, SC, and HN stations 

Visibility plots (Figure 2,3, and 4) and closure plots (Figure 
5,6, and 7) show that the SC and HN stations are calibrated, 
though not perfect. The visibility plots of HN-NL, and HN- 
SC shows quite a good match between the real amplitudes & 
phases and the calculated ones from the obtained map. The 
visibility plot of BR-SC shows good match in phase varia- 
tion, that of FD-HN shows good match in amplitude vari- 



ation. Presumably, the HN and SC stations are calibrated 
quite well. At least, the calibrations to the HN and SC sta- 
tions remove largely the systematic errors. 
About the MK station, the calibration solutions were ob- 
tained for the last one hour. Because the duration is short, 
the match in visibility is difficult to confirm. But, some of 
the triangles including MK station show nearly constant 
and/or very small variations in closure phase. These features 
presumably mean that the MK station got fringes between 
other stations though the signals arc weak. 
In APPENDIX D, we show several images by omitting the 
visibility data of the baselines to the MK, SC, HN stations. 
These images suggest that the calibration errors and the 
low SNR of the data from three stations have little influ- 
ence upon the image quality. 



2.5 Obtained image from the whole time data 

integration 

We made an image from the whole time integration using 
the task IMAGR in AIPS with loop gain parameters (GAIN) 
= 0.005, limit of subtracting iterations (NITER) = 20000, 
maximum residual flux density level (FLUX) = 0.005 Jy, 
and u-v Gaussian tapering (UVTAPER) = 6000 km. Here 
we used the wider tapering in u-v with expectations of being 
able to use lower SNR data points by the CALIB solutions. 
The mapping area is 3 mas (e — w) x 6 mas {n — s) and the 
boxing area is 2.13 mas (e — w) x 4.17 mas (n — s) at the 
center. Because the sizes of SgrA" at 4.'-)GHz arc w(>ll known, 
we limited the peak search to the central area to avoid se- 
lecting side lobe area. The grid numbers are 1024 (e — w) x 
2048 (n — s), namely the grid spacing is nearly 2.93 fias. By 
using a smaller grid, we intended to select the peak positions 
as accurately as possible. We performed the imaging with a 
restoring beam of 0.40 mas x 0.15 mas, PA = 0° 
The obtained image from the whole time integration sug- 
gests that the calibration of the visibilities arc success- 
fully performed. From direct fltting of an elliptical Gaussian 
shape to the image and the deconvolution of the restoring 
beam shape, we got a result of the shape and size. The ma- 
jor axis is 0.759t^;!555? mas, the minor axis is 0.323 ± 0.007 
mas, and the position angle is PA = 85.1 ± 1°. This estima- 
tion is consistent to the previous measurement of the same 
data using closure amplitude by Shen et al. (2005) that show 
the major axis is 0.722 ± 0.002 mas, and the minor axis is 
0.395t;!;o2o ™a,s and the PA is 80.4 ± 0.8 °. The deviations 
from the Shen's measurement are +0.037 mas in major axis, 
—0.072 mas in minor ax;is and 4.6° in PA. Miyoshi et al. (sub- 
mitted to PASJ) show the details of the results from other 
VLBA data on SgrA* at 43GHz with the same reduction 
manner. The consistent Gaussian shape and size to those 
found Shen et al. (2005) indicates that we got better cali- 
bration solutions from the CALIB and that we successfully 
obtained images of SgrA* with better spatial resolution than 
before. 

3 ANALYZING METHOD AND RESULTS 

We investigated the spatial distributions of oscillations in 
flux density with two independent methods. The first one is 
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Figure 5. The time variations of closure phases in respective triangles (1/3). 
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Figure 7. The time variations of closure phases in respective triangles (3/3). 




Figure 8. The top right panel (a) shows the Sgr A* image from the whole observing time. The contour levels are every 10 % of the 
peak brightness. We defined 209 small rectangles of 0.1 mas in east-west and 0.15 mas in north-south around the map center, distributed 
as matrix 19 in east-west by 11 in north-south. We measured the time variations of fiux densities of the rectangles. Panels (b), (c), and 
(d) show the flux density variations of the central three rectangles, 0.1 mas west, center, and 0.1 mas east along the central east-west 
line respectively. Panel (e) shows the time variations of the total of 209 rectangles. The rms noise levels, measured at the outer region of 
the maps whore no real emission should come are also shown as black points in Panel (d). The intensity integrated over 209 rectangles 
is about 1.3 Jy, which is the typical flux density of Sgr A* at 43 GHz, and those of the central 3 rectangles are 7 times higher than that 
of the rms noise level. Panel (f) shows the power spectra of periodicity calculated from the time variations using MEM. The x-axis is 
the period (minute) and the y-axis is the amplitudes (arbitrary unit). The blue line is that of the 0.1 mas-east rectangle, the green line 
is that of the central rectangle, and the red line is that of the 0.1 mas- west rectangle. Panel (g) shows the spectra of the flux density 
variation of all 209 rectangles. The black line in Panel (f) means the power spectra of the noise variations shown in Panel (d). The peaks 
of noise spectra appear at different periods from those of peaks of time variations in the central three rectangles. The bottom panel 
(h) shows the spreading size, one standard deviation of the distribution of each period of the area of the 209 rectangles. The black line 
shows the absolute size, while the blue line shows the component of the east-west direction. The red line shows that of the south-north 
direction. The five arrows show the positions of P = 12.9, 16.8, 22.2, 31.4, &;56.4 min from right to left respectively. At these periods, 
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Figure 9. The slit-modulation-imaging (SMI) maps with the 4 periods P = 16.8(a), 22.2(6), 31.4(c), and 56.4 min (d). The 8 maps 
show the images in phase frame series with respective periods. The white filled circles mark the intensity peak position in each map. 
In P = 31.4 & 56.4 min, second peaks are also marked with circles. The contour levels are every 10 % of the peak brightness of the 
respective maps. 



the comparisons of power spectra of time variations of inten- 
sities in 209 small rectangles areas of the image (Figure 8 a). 
We produced 360 snapshot maps from 5-minute integrations 
with a starting time increment of 30 seconds. Each map was 
produced by the same parameters as explained in the previ- 
ous section except the maximum residual flux density level 
(FLUX) =10mJy. 



We then measured flux density variation in each small 
rectangle. Their time variations of the flux densities show 
small up-and-downs similar to, but not the same as, those 
of other rectangles (Figure 8 b, c, d, and e). We then 
calculated their power spectra using a maximum entropy 
method (MEM, Burg 1967). In Figure 8 (f) are shown the 
power spectra of the central 3 rectangles, whose SNRs are 
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more than 7, while in Figure 8 (g) are the spectra of time- 
variations of all 209 rectangles shown in Figure 8 (e). 

Figure 8 (h) shows the spreading size cr{P), namely the 
one standard deviation of the spatial distribution of the am- 
plitude of each period. The definition is 

209 

a(P) = sgrt[^(4 + yt) ■ MP)/S{P)] (1) 

i=l 

The component in x direction is 

209 

(T.(P) = sqrtl^xf ■ Ai{P)/S{P)] (2) 

1=1 

The component in y direction is 

209 

ayiP) = sqrt[J2yi ■ MP)/S{P)] (3) 

i = l 

,where S{P) = J2Z\ MP), the sum of the amplitudes 
in the 209 rectangles. The i means i*^ rectangle, P is the 
period, Ai{P) is the amplitude of the period in each rect- 
angle, {xi, Ui) is the central position of the i*'' rectangle. 
If all of the power of the period concentrate at the central 
rectangle i = 0, the spreading sizes cr{P), ax{P), & crx{P) 
are zero. If the power distributes uniformly in all the rectan- 
gles, the spreading sizes are 0.362 mas in o-(P), 0.274 mas in 
ax{P), and 0.237 mas in ax{P)- In the case that the spread- 
ing sizes are larger than these values, the power of the period 
mainly comes from outer area, where the flux density and 
the SNRs are lower. Such periods with larger spreading sizes 
presumably originate from noise. If the spreading sizes of the 
periods are smaller than these values, the power of the pe- 
riods comes from the central part of the image, where the 
flux density and the SNR is higher. A real periodicity, if de- 
tected, because it should come from the central part of the 
image, the spreading size should be smaller than that of the 
uniform case. 

Small spreading sizes, or the concentrations of the 
distributions at the center, are seen around P = 
12.9, 17.2, 30.1, and 55.8 min, where the amplitude shows a 
maximum (Figure 8 g). Loose concentrations appear around 
P = 128.4, and 268 min though they are under-sampled. 

The first method tends to be affected by differences of u- 

V sampling between snap shot maps. We hence devised and 
used a new method, slit-modulation-imaging (SMI) which 
is a way of examining the existence of the periodic struc- 
tural change with an assuming period and almost free from 
u-v sampling differences (Miyoshi, 2008). SMI maps are ob- 
tained as follows: wo divided the whole observation time into 
several sections with the trial period P min interval and then 
each section was sub-divided into 8 phase-segments. The 
first-phase SMI map was produced with visibilities from the 
first phase-segments of all the sections. The second-phase 
SMI map was produced with visibilities from the second 
phase-segments of all the sections. Namely, the n-th phase 
SMI map is produced from the visibilities of the n-th phase- 
segments of all the sections. Thus, the SMI maps have u- 

V coverage very similar to each other, and are almost 
free from the effect of differing u-v coverages (Figure 10). 
Characteristic structure changes with the period will be en- 
hanced and emerge onto the respective SMI maps if it ex- 
ists. (Miyoshi 2008) 



Period 


Center 


Position 


Diameter 


PA 


Inclination 




RA 


DEC 








(min) 


(mas) 


(mas) 


(mas) 


(°) 


(°) 


16.8 


-0.02 


0.00 


0.177 


81.8 


87.7 


22.2 


-0.03 


-0.01 


0.149 


81.8 


87.0 



Table 1. Fitted parameters of circles to the displacements of the 
briglitest peaks in the P = 16.8 & 22.2 min 



We tested around these confined periods of P = 
12.9, 17.2, 30.1, and 55.8 min, adding with the P = 22.2 min 
detected at X ray flares (Belanger et al. 2006, Hamaus et al. 
2009) (Figure 8 (f) also shows a small peak at the period). 
Except for P = 12.9 min, we found conspicuous periodic 
change patterns from P = 16.8,22.2,31.4, and 56.4 min, 
with estimated errors ±1.4, ±1.4, ±1.5, ±6 min respectively. 
Because of the findings of simple change patterns, we regard 
the 4 periods as confirmed detection. They arc very roughly 
commensurable in a 3 : 4 : 6 : 10 ratio. At first glance these 
SMI results seem in stark contrast with the previous VLBA 
observations of SgrA*at 43GHz. For quantitative evaluation 
of the difference, we performed an elliptical Gaussian fitting 
directly to the obtained SMI maps, and the results of decon- 
volution of the restoring beam show that these SMI maps, if 
wc regard them as elliptical Gaussian shapes, are consistent 
to the previous size & shape estimation of SgrA* from clo- 
sure analysis (Figure 11). The differences are mostly within 
error bars of Bower's estimations (Bower et al. 2004) but 
out of the small error bars of Shon's estimations (Shen ct al. 
2005) 

In the case of P = 16.8 min, the peak is in the west of 
the image at the 1"* and the 2"'* phase frames, but shifts 
eastward, stays in the cast at the 3'^'', 4*'' and 5*'' phase 
frames, and after that, moves westward again (Figure 9 a). 
In the case of P = 22.2 min, at the 1"' phase frame the 
peak is in the east, moves towards the west at the 2"'* 
phase frame, shifts further west at the 4*^ and 5*^ phase 
frames, and then moves back towards the east from the 
Qth gth pjja,se frames (Figure 9 b). Namely, in the 
P = 16.8 and 22.2 min structure changes, the brightest 
position oscillates mostly in the east-west direction with 
the respective periods, which looks like a one-arm structure 
(m = 1) motion viewed edge-on. 

Assuming the motions to be projected circular ones, we 
found the P = 16.8 min orbit is 0.18 mas while that of the 
P = 22.2 min is 0.15 mas in diameter, both are tilted with 
PA = 82° nearly the same position angle as the major axis 
of the average SgrA*imagc at 43 GHz (Figure 12). The 
viewing angles are almost edge-on, 87.0° (P = 16.8min) 
and 87.7° (P = 22.2 min) respectively (Table 1). 

In the P — 31.4 min change, one or two peaks appear 
in each phase frames (Figure 9 c). The numbers of peaks 
arc 2— >1— >2— >2^1— >2— >1— >1 with phase frame 
change. At first glance this seems random, but a relation 
exists. The sum of peak numbers in the two SMI maps just 
half a period away from each other is always three, which 
can be explained by the existence of a three-arm pattern (m 
= 3) with evenly separation angle rotating with the period. 
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(a)P = 16.8min (b) P = 22.2 min (c) P = 31.4 min (d) P = 56.4 min 
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Figure 10. The u-v coverage of the respective SMI maps. 
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Frame Frame Frame Frame 

Figure 11. Gaussian size estimations of the each SMI maps. The dashed Une shows the parameter obtained from the whole time 
integration image. The dot line shows the parameters estimated from closure analysis (Shen et al. 2005). 



which we view edge-on. 

In the P = 56.4 min change, one peak appears during 

a half period (the 6^^ to S***, and 1"* phase frames) and two 
peaks appear during the last half period from the 2"'' to the 
S*** phase frame (Figure 9 d). This is due to the amplitude of 
variation at the central 50 /uas being twice as large as those 
of the surroundings. One peak appears when the central part 
becomes bright, and two peaks, one on the east and one on 
the west, appear when the central part darkens. Unlike other 
periods, the P = 56.4 min component seems concentrated 
and excited at the radius close to the black hole where the 
relativistic effect are not negligible, and hence it can be a 
key to distinguishing the origin of QPOs. 

The detected pattern changes are not from noise ef- 
fects: we simulated changes using artificial visibilities with 
noise levels showing the comparable amplitude fluctuations 
with those of real visibility data. SMI images from the arti- 
ficial visibilities show larger variations in not east-west but 
south-north direction that the spatial resolution is worse. 



while those from the observations show the opposite ten- 
dency. From the artificial visibilities, changes appeared in 
the image but none of the periodic regular patterns could 
be found. Also the artificial visibilities showed no coincident 
occurrence between spreading size minima and amplitude 
peaks in the spectra of the whole 209 area. Further, the ob- 
served periodic changes of real data basically did not disap- 
pear in SMI- maps with shorter total integration times (195 
min) while different changes appeared from the artificial vis- 
ibilities. 



4 DISCUSSION AND CONCLUSION 

The radio QPOs reported here have a delay of 1.5 days from 
the millimeter wave flare, and the distributions extended to 
100 -Rsin apparent diameter while the QPOs in NIR and 
X ray are just at the flaring time and estimated to occur 
around a marginally stable orbit (MSO) (Genzel et al. 2003; 
Aschenbach et al. 2004; Eckart et al. 2006 and Belanger et al. 
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Figure 12. The offset of brightness pealc positions in the frame series of SMI images: those of P= 16.8 min in the left panels and those 
of P= 22.2 min in the right panels. The top panels show the positional shifts in R.A. while the middle panels in declination. The curve 
in each panel shows the sinusoidal function obtained from fitting to the data. The bottom panels show the displacements of the peaJjs in 
the sky. The ellipses are formed from the obtained fitted sine curves. 



2006). The flaring events at MSO would vibrate the whole 
system to excite fundamental disk oscillation modes, which 
would last for several days at least. 

There arc features that cannot be explained as Kcplc- 
rian motions directly governed by centripetal force, but that 
are easily interpretable as disk oscillation features (Kato et 
al. 2008). First, the counter-rotations occurred simultane- 
ously between the P — 16.8 min and P = 22.2 min periods, 
which is difficult to be explained by a single body motion 
like hot spot models. Second, the orbit of the shorter period 
P = 16.8 min is larger than that of the P = 22.2mm, which 
cannot bo realized by a Keplerian motion. 

If QPOs originate in a strong gravity field where the 
relativistic effect plays an important role, the periods of 
QPOs should depend on the mass and the spin of a massive 
black hole. Recent theories of disk seismology (Kato et 
al. 2008; Abramowicz et al. 2001) predict that pealc fre- 
quencies of QPOs can be scaled by a mass of central black 



holes as an analogy to QPOs in black hole x-ray binaries 
(BXB) (Remillard, & McClintock 2006). For example, in 
GRO J1655-40, a peak frequency of high frequency QPOs is 
about 3 X 10^(6.0 - 6.6Mq/Mbh) Hz (where Mq is a solar 
mass), with the result that a corresponding peak frequency 
using tlic mass of Sgr A* derived from the orbital motions of 
surrounding stars (3.6 ± 0.3 x IO^Mq) (Eisenhauer et al. 
2005) is about 5.1 x 10"" Hz (P = 32min), which is one of 
our findings. Detailed analysis with the obtained four QPO 
periods and wave-warp resonant oscillation model predicts 
the spin of SgrA*to be 0.44 ± 0.08 and the black hole mass 
to be (4.2 ± 0.4) x 10® Mq. (Kato et al. submitted to MN) 

The apparent angular size of 0.1 mas corresponds to 
10.4_R,5 (l_Rsis Schwarzschild radius) assuming the Galactic 
center distance to be 7.6 kpc. (Eisenhauer et al. 2005), 
which scale give us that the orbital velocities of bright 
peak positions are vp=i6.8min = 2.1 c at 8.8 Rs and 
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vp=22.2min ~ 1-4 c at 7.4 Rs, both are superluminal, 
but can be explained as follows. It is well-known that 
the VLBI images of SgrA* are obscured and broadened 
by the scatterings of surrounding plasma (Lo et al. 1998; 
Doeleman et al. 2001). If we consider the broadening 
ratio of 2.6 ~ 3.0 at 43 GHz by scattering (Shen et 
al. 2005; Bower et al. 2004), and the magnifying ratio 
~ 1.23 by self-gravitational lensing effect at 3 -Rs, the 
apparent 0.1 mas observed at 43 GHz is equal to 2.8 to 
3.3 -Rsin the intrinsic image of SgrA*. Taking into these 
effects, we find that Wp^jg = 0.43 c at 2.7 Rs and 

i''p=22.2min ~ 0.68 c at 2.3 Rs , which velocities are 
comparable to the Keplerian velocities at r = 2 — ARs from 
a black hole. 

Presumably the emission region of our observed radio 
QPOs is different from those of other wavelength, even 
though the origin of oscillations is common. Kato et al. 
(2009) show that synthetic images of radio, NIR, and X-ray 
are quantitatively different, by using multi-wavelength 
radiative transfer calculations including compton scattering 
process in three-dimensional magnetized hot accretion 
flows. Without synchrotron self-compton (SSC) process, 
Takahashi et al. (2009) show that temporal fluctuation of 
flux density at shorter wavelength, where the disk becomes 
optically thin, is insensitive to the non-axisymmetric disk 
structure. A possible emission mechanism of NIR and X-ray 
wavelengths in order for having common oscillations at 
radio wavelength is the SSC, which is dominant when mass 
accretion rate becomes large [see Fig. 4 (a) and (d) in Kato 
et al. 2009]. 

The amplitude of the QPOs were 22 mjy {P = 
16.8min), 24mJy (P = 22.2 min), 25mJy (P = 31.4min), 
and 20 mJy (P = 56.4 min) respectively, the total reached 90 
mJy, nearly 10 % of the constant flux density component of 
the image {— 980 mJy). The relative intensities of their pe- 
riods are about 2 or 3% level of the total flux density, which 
is comparable ratio to those of X ray QPOs (5 — 10% level in 
case of X ray Low Frequency QPO, 1% level in case of X ray 
High Frequency QPO). Such structural changes will cause 
fluctuations in fringe phase, rate and amplitude. Existence 
of oscillations with periods shorter than the whole observ- 
ing time will be one of the reasons why the VLBI imaging 
of SgrA* is difficult so far. However at the same time, the 
detection of these periodic pattern changes means the ob- 
scured radio image of SgrA* still bears the imprint of the 
intrinsic figure. The VLBA observations of SgrA*, together 
with oscillation analysis, will provide us a unique chance to 
investigate the strong gravity field around the massive black 
hole SgrA* with order of Pgresolutions. 
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APPENDIX A: THE SMI RESULTS AROUND P 
= 12.9 MIN 

We here show the none detection cases of conspicuous pat- 
tern change with the trial periods using the SMI method. 
As shown in Figure 8 (f), (g), there is a small peak at the 
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P=12.2min P = 12.4 min P=12.6min P = 12.8 min P=12.9min P = 13.0 min P=13.1min P = 13.2 min P = 13.4 min 




Figure Al. SMI search around P = 12.9 min. The sUt-modulation-imaging (SMI) maps with the 9 periods P = 
12.2, 12.4, 12.6, 12.8, 12.9, 13.0, 13.1, 13.2, an(il3.4min. The 8 maps show the images in phase frame series with respective periods. The 
cross mark shows the position of peaJc intensity in each map. The contour levels are every 10% of the pealc brightness of the respective 
maps. Prom these trial periods we detected no conspicuous pattern change with the period. 



period P — 12.9 min in the power spectra. Also the Figure 
8 (h) shows that the spreading size of the P = 12.9 min is 
relatively small. Using the SMI method we searched around 
the period for any conspicuous pattern changes. One-peak 
or two-peaks structures appeared in the phase frames of SMI 
maps of the respective trial periods. Unlike the cases of the 
P = 16.8, 22.2, 31.4, and 56.4 min, the change patterns in 
the phase frames seem at random. We therefore could not 
confirm the QPO detection around P = 12.9 min. 



APPENDIX B: 
ARRAY 



SPECTRA FROM SMALLER 



The QPO spectra can be detected from images with 
limited u-v coverage. In order to indicate that the obtained 
QPO spectra are not due to insufficient calibrations of 
the stations of MK, SC, and HN, we produced series of 5 
minutes integration snap shot maps without the visibilities 
of the longer baselines, and performed the same analysis of 
the time variations of flux densities in several regions of the 
mapping area. 

We omitted the following baselines from the original 
calibrated data. 

• All the baselines to the MK station. 

• The 7 baselines to the SC stations, namely the baselines 
of SC-BR, SC-FD, SC-KP,SC-LA, SC-NL, SC-OV, and SC- 
PT. We keep the data of the SC-HN baseline because the 
visibilities of SC-HN show quite high SNR. 



• The 3 baselines to the NL station, namely NL-BR, NL- 
OV, and NL-SC baselines. 

It is already noted in the section 2 that the original cal- 
ibrated data do not contain the baselines of MK-SC, and 
MK-HN. Figure Bl (b) shows the u-v coverage of the limited 
data set. For comparison, we show the whole u-v coverage 
of the original calibrated data in Figure Bl (a), which is the 
same figure as Figure 1. 

Figure B2 (a) is the same image as shown in Figure 8 (a) 
with a different contour expression, while Figure B2 (b) 
shows the image of the limited data set from the whole time 
integration. We here used the common restoring beam of 
OAOmas x O.lbmas, PA — 0°. Other parameters of the 
imaging are the same as written in the section 2-5. Be- 
cause of the poor u-v coverage of the limited data, the re- 
sultant image is a simple obuje only representing the size 
and the elongation direction (east-west) of the SgrA'shape 
at 43GHz. 

From the series of snap shot maps, we measured the 
flux densities in the 8 regions defined as shown in Figure 
B3, and analyzed the time variations using MEM and found 
the power spectra of periodicity shown in Figure B4. There 
are two facts to be noted. First, both of the spectra from the 
whole u-v data (black lines) and the limited u-v data (red 
lines) show quite similar curve to each other in respective 
regions. Second, peaks appear around P = 17, 22, 30, 55 
min in the spectra of the limited u-v data. These two facts 
indicate at least that the detected QPOs are not due to 
insufficient calibrations of the stations of MK, SC, and HN. 
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Figure Bl. The u-v coverages of the whole caUbrated data (a) and the hmited data (b). The Figure B 2 (a) is the same plot as Figure 
1. 
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Figure B2. maps 



APPENDIX C: CALCULATED CLOSURE 
PHASES FROM THE SMI MAPS (P =16.8 MIN) 

In Figure C 1 and 2, wc show the closure phases calculated 
from the 8 maps from the SMI method at P= 16.8 min. 
The SMI maps show consistent closure phases to the ob- 
served ones. First of all, the small triangles containing FD, 
LA, PT, KP, and OV stations shows nearly zero, constant 
closure phases both in observed closures and the calculated 
ones. Second, the observed closure phases of triangles con- 
taining one station from BR or NL and two stations from the 
inner five stations FD, LA, PT, KP, and OV show deviations 



of a few tens of degrees from zero. The calculated closure 
phases from the SMI maps show different values from the 
observed ones, however, the distributions of them cover the 
regions of variations of the observed closure phases. Third, 
the larger the triangles become, the observed closure phases 
show larger deviations up to ±180°. The calculated closure 
phases from the SMI maps also show the same large devia- 
tions. From the point of closure phases, the SMI maps show 
consistent behavior to the observed closure phases. 
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Figure B3. The definition of the 8 regions. The size of region (4) is comparable to so called a scattering size of SgrA'at 43GHz. Assuming 
the distance to the Galactic center to be 8 kpc and the mass of SgrA*to be 4 X Msun, the apparent 3 mas corresponds to be 24 AU, 
or 305 Schwarzschild radii. 
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Figure B4. The power spectra of the periodicities in the 8 regions. Black line shows the spectrum obtained from the whole calibrated 
data, red line shows the spectrum obtained from the limited calibrated data. 
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Figure C2. The calculated closure phases from the eight SMI maps (P =16.8 min) (2/2) 
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APPENDIX D: IMAGES WITH SEVERAL 
LIMITED VISIBILITY SETS 

Figure D shows sc-voral images from the whole time integra- 
tion with several limited visibility sets. Here we performed 
imagings with the task IMAGR with the same parameters 
as described in the section 2-5. Visibility data from the base- 
lines to MK seem to do almost nothing to the image quality. 
The 4 images without the baselines to MK, (b), (d), and (e) 
in Figure D, are very similar to each other and also to the 
image from all data, (a) in Figure D. While the existences of 
the baselines to HN, and SC in data play fairly an important 
role to the imaging quality. The exception of one from the 
two stations gives little influence to the image quality ((c), 
(d), (e), (f), (g) in Figure D). The exception of both of the 
HN, SC stations gives quite a degree of influence in image 
quality ((g), and (h) in Figure D). 
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Figure Dl. Several images with limited visibility data, (a) the image using all the calibrated data, the same one as shown in Figure 8 
(a), (b) an image using all the calibrated data except baselines to MK station, (c) an image using all the calibrated data except baselines 
to HN station, (d) an image using all the calibrated data except baselines to SC station, (e) an image using all the calibrated data except 
baselines to MK and SC stations, (f) an image using all the calibrated data except baselines to MK and HN stations, (g) an image using 
all the calibrated data except baselines to SC and HN stations, (h) an image using all the calibrated data except baselines to MK, SC, 
and HN stations. 



